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Beschreibung 

Rechnergesttltztes Verfahren zur parallelen Berechnuiig des Ar- 
beitspunktes elektrischer Schaltungen 

5 

Die rechnergestUtzte Simulation elektrischer Schaltungen hat 
bei der Entwicklung sehr groBer Schaltungen, also Schaltungen 
mit einer sehr groBen Anzahl von Elementen, eine immer grofle- 
re Bedeutung erlangt. Insbesondere bei der Entwicklung von 
10 Computerchips mit einer Vielzahl, beispielsweise mehrere hun- 
derttausend Transistoren, erweist sich eine serielle Bearbei- 
tung fiir die Ermittlung der SchaltungsgroJien durch einen 
Rechner wegen des zu groiien Zeitauf wands als unbrauchbar. 

15 In der WO 98/24039 wird deshalb vorgeschlagen, eine groBe 

Schaltung zu partitionieren und die Partitionen jeweils von 
verschiedenen Rechnern bearbeiten zu lassen. 

Oblicherweise werden bei der Berechnung zuerst der Arbeits- 
20 punkt, das heiBt die Potenziale aller Knoten als Grundlage 
ftlr weitere Analysen, wie beispielsweise Transienten- oder 
Wechselstromanalysen, ermittelt. 

Zur parallelen Berechnung wird in U. Wever, Q. Zheng et al.: 
25 "Domain Decomposition Methods for Circuit Simulation" 

(Proceedings of the 8^^ Workshop on Parallel and Distributed 
Simulation, PADS '94 Edinburgh, Scotland, UK, S. 183-186, Ju- 
li 1994) und in U. Wever, Q. Zheng: "Parallel Transient Ana- 
lysis for Circuit Simulation" (Proceedings of the 29^^ Annual 
30 Hawaii International Conference on System Sciences, S. 442- 
447, 1996) eine Implementierung des Newton-Verf ahrens vorge- 
schlagen. Nachteiligerweise kann hier, aufgrund schlechter 
Konvergenzeigenschaf ten, nur bei Vorliegen hinreichend guter 
Schatzungen des Arbeitspunktes Konvergenz erzielt werden, 
35 Derartige gute Schatzungen sind in der Regel jedoch bei gro- 
Ben Schaltungen schwer oder nicht zu bewerkstelligen. 
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Der vorliegenden Erfindung liegt deshalb die Aufgabe zugrun- 
de, ein rechnergesttitztes Verfahren zur parallelen Berechnung 
des Arbeitspunktes elektrischer Schaltungen zu schaffen, das 
eine einfache, sichere und schnelle Berechnung des Arbeits- 
5 punktes der Schaltiing gewahrleistet . 

Diese Aufgabe wird erf indungsgemafi durch das Verfahren gemalJ 
Patentanspruch 1 gelost, 

10 Bei dem Verfahren wird zur parallelen Berechnung des Arbeits- 
punktes das Auf ladeverf ahren verwendet, welches sehr gute 
Konvergenzeigenschaf ten auf weist . 

Durch das Auf ladeverf ahren werden, wie beispielsweise in H. 

15 Spire: "Simulation integrierter Schaltungen" (2. Auflage, R. 
Oldenbourg Verlag Mtlnchen, Wien 1990) beschrieben, die in ei- 
ner elektrischen Schaltung vorhandenen dynamischen Elemente, 
wie Kapazitaten und/oder Induktivitaten, dazu verwendet, ura 
ilber eine Pseudo-Transientenanalyse den Arbeitspunkt zu be- 

20 rechnen. Die dynamischen Elemente werden dabei schrittweise^ 
von einem Wert Null ausgehend, auf einen im Prinzip beliebig 
hohen Wert, insbesondere auf den Wert "1" hochgef ahren. 

In einer Ausftihrungsf orm nach der Erfindung konnen zusatzlich 
25 an geeigneten Stellen, beispielsweise an oder zwischen mehre- 
ren Knoten dynamische Elemente, wie beispielsweise Kapazita- 
ten Oder Induktivitaten, mit vorbestimmten bekannten Werten 
zusatzlich eingebaut werden. Hierdurch konnen Schwierigkeiten 
bei der Berechnung aufgrund einer zu geringen Anzahl von dy- 
30 namischen Elementen vermieden werden. 

Hierzu konnen die zusatzlich in die Schaltung eingebrachten 
dynamischen Elemente von einem 7\nf angswert, beispielsweise 
gegen Unendlich, fur den die Arbeitspunktberechnung trivial 
35 ist, schrittweise auf Null verringert werden, so dass wieder- 
um die ursprtingliche Schaltung simuliert wird. 
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In eiher bevorzugten Ausftihrungsf orm der Erfindung wird we- 
nigstens ein Knotenpunkt/ also ein Verbindungspunkt von we- 
nigstens zwei Strompfaden^ der Schaltung mittels einer Kapa- 
zitat mit einem vorbestimmten Wert mit einem vorbestimmten 
Potenzial verbiinden. 

In weiterer Ausgestaltiing der Erfindung kann aber auch an je- 
den Knoten xiber alle Partitionen eine Kapazitat angeschlossen 
werden, deren zweiter Anschluss jeweils an einem vorbestimm- 
ten Potenzial, beispielsweise an Masse liegt. Diese Vorge- 
hensweise hat den Vorteil, dass der Arbeitspunkt, also die 
jeweiligen Potenziale der einzelnen Knotenpunkte, fur einen 
gegen Unendlich gehenden Anfangswert fur die Kapazitaten, 
aufgrund einer Losung der dann expliziten Gleichung, fiir die 
Schaltung trivial zu berechnen ist. Ober eine schrittweise 
Veranderung des Wertes ftlr die Kapazitaten kann dann durch 
geeignete Neuwahl des Wertes fUr die Kapazitaten die Schal- 
tungssimulation verandert werden, bis sich eine Berechnung 
des Arbeitspunktes der Schaltung ftlr einen Wert der Kapazita- 
ten gegen Null oder nahezu gegen Null ergibt. 

Durch diese Implement ierung des Auf ladeverf ahrens auch far 
parallele Berechnung einer elektrischen Schaltung kann vor- 
teilhafterweise eine sehr grofie Schaltung mit einer Vielzahl 
von Transistoren selbst bei einer geringeren Zahl von dynami- 
schen Elementen auf schnelle und einfache Weise berechnet 
werden. Ftlr die jeweilige Neuf estlegung des Wertes fiir die 
Kapazitat mit dem Ziel, diesen gegen Null gehen zu lassen, 
sind verschiedene Vorgehensweisen denkbar, wobei als Ent- 
scheidungskriterium der Schwierigkeitsgrad der Berechnung des 
Arbeitspunktes des jeweiligen vorhergehenden Schrittes in 
Frage kommt. 

Der Arbeitspunkt, der sich durch Losung einer nicht-linearen 
Gleichung ftlr jeweils einen bestimmten Wert der Kapazitaten 
ergibt, kann beispielsweise iterativ mit dem Newton-Verf ahren 
gelost werden. Ftlr die Wahl des nachsten Wertes fur "C" kon- 
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nen dann die Anzahl der zur LSsung notwendigen Iterations- 
schritte verwendet werdeii/ bis ein vorbestimmter Wert ftlr "C" 
unterschritten wird. 

Weitere vorteilhafte Ausgestaltungen der Erfindung ergeben 
sich aus den Unteranspriichen. 

Die Erfindung wird nachfolgend anhand eines in der Zeichnung 
dargestellten AusfUhrungsbeispiels erlauter. In der Zeichnung 
zeigt die einzige Figur ein Ablauf diagramm, in dem die ein- 
zelnen Verf ahrensschritte des Verfahrens dargestellt sind. 

Nach der Zeichnung liegt als Ausgangspunkt 1 eine Simulation 
einer elektrischen Schaltung, beispielsweise in der Schal- 
tungsbeschreibungssprache SPICE vor. Diese wird in einem er- 
sten Verarbeitungsprozess 2, wie beispielsweise in WO 
98/24039 als "Clustering-Verf ahren" beschrieben, partition 
niert/ so dass sich einzelne Partitionen beziehungsweise Tel- 
le der Schaltungen^ die moglichst mit gleichem Schwierig- 
keitsgrad berechenbar sind, ergeben. 

Selbstverstandlich sind auch andere anwendbare Partitionie- 
rungsvarianten, wie beispielsweise die in N. Frohlich, B. 
Riess, U. Wever, Q. Zheng: "A new approach for parallel simu- 
lation of VLSI circuits on a transistor level" (IEEE Transac- 
tion on Circuits and Systems - I: Fundamental Theory and 
Applications, Vol. 45, No. 6, June 1998, Seiten 601 bis 613) 
beschriebene "Ratio-Cut -Methode" Oder selbst eine willkurli- 
che Aufteilung "per Hand" denkbar. 

In einem weiteren Arbeitsschritt 3 wird nun jedem Knotenpunkt 
der Schaltung, also einem Verbindungspunkt von wenigstens 
zwei Leitern beziehungsweise Strompfaden eine geerdete Kapa- 
zitat hinzugeftigt. Selbstverstandlich ist es auch denkbar, 
eine Kapazitat, deren zweiter Anschluss an einem vorbestimm- 
ten Potenzial liegt, hinzuzuftigen, wobei sich sowohl Pdten- 
ziale als auch die Werte fur die Kapazitaten jedes Knoten un- 
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terscheiden k5nnen. Aus Effizienz-, StcJDilitats- und berech- 
nungs6konomischen Grtlnden liegt jede Kapazitat mit ihrem 
zweiten Anschluss an dem gleichen Potenzial/ beispielsweise 
Masse, wobei zudem aus denselben Grtinden ftir alle Kapazitaten 
5 in Schritt 4 ein identischer hoher Wert CO gewahlt wird. 

Fur diesen Wert CO wird nun in einem weiteren Schritt 6 fiir 
jede Partition beziehungsweise Teilschaltung der Arbeitspunkt 
berechnet/ wobei die erforderlichen Kopplungswerte, also die 
10 Werte fiir die Kopplungspunkte beziehungsweise Schnittstellen 
benachbarter Partitionen, ftir die Berechnung der Arbeitspunk- 
te benachbarter Partitionen ausgetauscht und einberechnet 
werden. 

15 Hierbei kann aus Ef f izienzgrtinden vorteilhaf terweise eine 

Partition, welche im Folgenden als "Master" bezeichnet wird, 
die Steuerung des Auf ladeverf ahrens tibernehmen, Der Master 
bestimmt dann den Anfangswert CO fur die Kapazitaten, wobei 
es auch denkbar ist, den Anfangswert extern, beispielsweise 

20 vom Benutzer, vorzugeben. Dieser Wert CO wird dann an alle 
Obrigen Partitionen, im Folgenden als Slaves bezeichnet, 
tibergeben, Nachfolgend wird sowohl in der Master-Partition 
als auch in alien Slave-Partitionen der Arbeitspunkt berech- 
net, wobei selbstverstandlich der voile Quellenvektor, der 

25 die in der Schaltung vorhandenen Energiequellen darstellt, an 
der Schaltung anliegt. 

Ausgehend von einem zu losenden Gesamtproblem 

X = (m,5,,...,5^) 

30 wobei m die Unbekannten des Masters, 
p die Anzahl der Partitionen, 
Si die Unbekannten der Partition i 

bezeichnet, fUhrt das Auf ladeverf ahren hierbei auf das Diffe- 
35 rent ialgleichungssys tern. 
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wobei 

D=diag(C, . . .C, 0, . . . , 0) und 
t die Zeit bezeichnet- 

Dieses Dif f erenzialgleichungssystem kann beispielsweise mit 
dem impliziten Euler-Verf ahren gel5st werden^ was auf die 
nicht-lineare Gleichung 



fuhrt. Dabei ist die Losung zum Zeitpunkt fur k = 0, 1^ 
usw. und h die Schrittweite t*^"^^ - t^. Diese nicht-lineare 
Gleichung kann beispielsweise iterativ mit deiti Newton- 
Verfahren 



fUr n = 0, 1, usw. gelost werden, 

Bei der parallelen Berechnung werden dann vom Master nur die 
Unbekannten "m" berechnet, wobei die Unbekannten "si" fUr den 
Master Festwerte darstellen. Damit eine Losung dieses berech- 
net werden kann^ werden von den Slaves bzw. Slave-Partitionen 
in jedem Iterationsschritt des Masters die korrigierten Werte 
berechnet und dem Master mitgeteilt . Hierzu muss von den Sla- 
ves das nicht-lineare Gleichungssystem 
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gel5st werden, wobei fi und Di die entsprechenden Gleichungen 
und Matrizen ftir den Slave "i" bezeichnen. Diese nicht- 
lineare Gleichung kann ebenfalls mit dem Newton-Verf ahren be- 
rechnet werden, wobei nicht notwendigerweise alle Iterations- 
schritte bis zur Konvergenz durchgeftihrt werden mtissen. 

Die Ergebnisse der Slaves werden daraufhin in das Gleichungs- 
system des Masters eingesetzt, woraufhin der Master den Aus- 
druck m^ll berechnen kann. 

In einem nachsten Schritt 7 wird ein neuer Wert "Cneu" fur "C" 
vom Master bestimmt, wobei auch die Slaves Vorschlage machen 
kSnnen. Um das Ziel, namlich einen Wert ftir "C" gegen Null 
beziehungsweise kleiner als einen vorbestimmten Wert e und 
damit eine Berechnung der ursprunglichen Schaltung zu errei- 
chen, wird die Wahl des neuen Wertes "Cneu" ftir "C" in Abhan- 
gigkeit von der Schwierigkeit der Berechnung des vorhergehen- 
den Schrittes mit dem Wert "Cait" ftir "C" getroffen. 

Hierbei sind verschiedene Vorgehensweisen denkbar, beispiels- 
weise die Wahl von "Cneu" in Abhangigkeit einer Analyse der 
Anzahl der Iterationsschritte, die der Master zur Berechnung 
des nicht-linearen Gleichungssystems ftir den vorhergehenden 
Wert "Cait" ftir "C" benotigt hat: 




72 n < ni 

C,^ fUr ni ^ n ^ n2 
fUr n2 > n 
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Wobei "n" die Anzahl der Iterationen des Masters und "xii"^ 
"n2" vom Benutzer vorgegebene Parameter sind. Anstelle der 



verstandlich auch andere Strategien zur Verkleinerung bezie- 
hungsweise VergroBerung von "Cait" inoglich. Zusatzlich kann 
auch die Anzahl der Iterationsschritte, die die Slaves bzw. 
Slave-Partitionen zur Losung ihres nicht linaren Gleichungs^ 
systems benotigt haben, beriicksichtigt werden, beispielsweise 
durch die Wahl von 



wobei ni, siavej die Anzahl der Iterationsschritte des Slaves j, 
wahrend der i-ten Iteration des Masters bezeichnet, Hierbei 
ist anzumerken, dass die Wahl eines Wertes far "C" nur einen 
Einfluss auf die Effizienz des Verfahrens, nicht jedoch auf 
den Arbeitspunkt selbst hat. 

Nach Erreichen eines Wertes von fUr "C" kleiner oder gleich 
einem vorgegebenen Wert e wird die Berechnung an der Verzwei- 
gung 5 abgebrochen, wobei in eineiti letzten Schritt der Wert 
fur C auf "0" gesetzt werden kann. Hierdurch steht als Ergeb- 
nis 8 der Arbeitspunkt der ursprunglichen Schaltung fest und 
kann uber Ausgabeeinheiten^ wie beispielsweise Bildschirm, 
Drucker oder ahnlichem, ausgegeben und/oder in einem Speicher 
als Grundlage ftir weitere Analysen der Schaltung zwischenge- 
speichert werden. 

Durch das erf indungsgemafie Verfahren kann ein Arbeitspunkt 
einer sehr groBe elektrische Schaltung vorteilhaf terweise von 
einer Vielzahl von Rechnern oder Prozessoren parallel berech- 
net werden, wobei der Nachteil bekannter paralleler Berech- 
nungsarten, namlich fehlende Konvergenz bei ungunstigen An- 
fangswerten, vermieden werden kann. Der Auf wand fur die Be- 
rechnung hinsichtlich Iterationsschritten und schrittweisem 



Halbierung beziehungsweise Verdopplung von "Cait" sind selbst- 
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Setzen eines Wertes ftlr "C" wird durch die parallels Berech- 
nung aufgriind der Verteilung auf mehrere Prozessoren in einem 
vertretbaren Rahmen gehalten. 

Der so berechnete Arbeitspunkt dient dann als Grundlage filr 
weitere Analyseii/ beispielseise der Wechselstromanalyse einer 
Schaltung. 
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Patentansprtiche 

1. Rechnergesttitztes Verfahren zur parallelen Berechnung des 
Arbeitspunktes elektrischer Schaltungen, 

5 - bei dem die Schaltung in einem ersten Schritt in mehre- 

re Partitionen partitioniert wird, 
dadurch gekennzeichnet, dass 
- zur parallelen Berechnung der einzelnen Partitionen das 
Auf laderverf ahren angewendet wird. 

10 

2. Rechnergesttitztes Verfahren nach Anspruch 1, 
dadurch gekennzeichnet, dass an we- 
nigsten einem Knotenpunkt der Schaltung ein dynamischess 
Element (C, L) vorgesehen wird. 

15 

3. Rechnergesttitztes Verfahren nach Anspruch 2, 
dadurch gekennzeichnet, dass an je- 
dem Knotenpunkt der Schaltung ein dynamisches Element (C, 
L) vorgesehen wird. 

20 

4. Rechnergesttitztes Verfahren nach Anspruch 2 oder 3, 
dadurch gekennzeichnet, dass jeder 
Knotenpunkt der Schaltung mittels jeweils einer Kapazitat 
mit jeweils einem vorbestimmten Wert mit jeweils einem 

25 Potenzial verbunden wird, so dass ein Arbeitspunkt der 

modif izierten Schaltung berechnet werden kann, 

5. Rechnergesttitztes Verfahren nach Anspruch 4, 
dadurch gekennzeichnet, dass 

30 an jedem Knoten einer Partition eine Kapazitat mit dem- 

selben Wert (CO) vorgesehen wird.. 

6. Rechnergesttitztes Verfahren nach Anspruch 4 oder 5, 
dadurch gekennzeichnet, dass 

35 jeder Knoten einer Partition mittels einer Kapazitat mit 

demselben Potenzial verbunden wird. 



•1 
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7. Rechnergesttltztes Verfahren nach Anspruch 4, 
dadurch gekennzeichnet/ dass 

an jedem Knoten aller Partition eine Kapazitat mit dem- 
selben Wert (CO) vorgesehen wird. 

5 

8. RechnergestUtztes Verfahren nach Anspruch 4 oder 7, 
dadurch gekennzeichnet, dass 
jeder Knoten aller Parti tionen mittels einer Kapazitat 
mit demselben Potenzial verbunden wird. 

10 

9. RechnergestUtztes Verfahren nach einem der Anspriiche 4 
bis 7, dadurch gekennzeichnet, dass 
das Potenzial auf Masse liegt. 

15 10. Rechnergesttltztes Verfahren nach einem der Anspruche 4 

bis 8, dadurch gekennzeichnet, dass 
- bei geeigneter schrittweiser Veranderung des Wertes (C) 
der Kapazitat jeweils der Arbeit spunkt der Schaltung 
berechnet wird und 
20 - dieser Schritt solange wiederholt wird bis die Werte 

der Kapazitaten nahezu Null sind. 

11. Computerprogrammprodukt, welches in einen Arbeitsspeicher 
einer Rechenanlage geladen werden kann, mit einem Softwa- 

25 recede fur das DurchfUhren des Verfahrens nach einem der 

vorhergehenden Ansprtlche, wenn das Computerprogrammpro- 
dukt auf einer Rechenanlage lauft. 

12. Datentrager mit einem Computerprogrammprodukt nach An- 
30 spruch 11. 
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Zuscunmenf assung 

Rechnergesttitztes Verfahren zur parallelen Berechnung des Ar- 
beitspunktes elektrischer Schaltungen 

5 

Die Erfindung betrifft ein rechnergesttitztes Verfahren zur 
parallelen Berechnung des Arbeitspunktes elektrischer Schal- 
tungen, bei dem die Schaltung in einem ersten Schritt in meh- 
rere Partitionen. partitioniert wird, bei dem zur parallelen 
10 Berechnung der einzelnen Partitionen das Auf laderverf ahren 
angewendet wird. 



15 



Hauptzeichnung ist die einzige Figur. 
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